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Abstract 

The helical Dirac fermions at the surface of topological insulators show a strong circular dichroism 
J> which has been explained as being due to either the initial-state spin angular momentum, the 

\Q initial-state orbital angular momentum, or the handedness of the experimental setup. All of these 



interpretations conflict with our data from Bi2Te3 which depend on the photon energy and show 

co 
o 

CO 



CO 

O several sign changes. Our one-step photoemission calculations coupled to ab initio theory confirm 



the sign change and assign the dichroism to a final-state effect. The spin polarization of the 

. — , photoelectrons, instead, remains a reliable probe for the spin in the initial state. 

>< 
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Topological insulators are characterized by an insulating bulk and a metallic surface. The 
gap in the bulk bands is due to a band inversion caused by strong spin-orbit interaction [IHZ] ■ 
Typical examples for which these gaps are rather large are the systems Bi 2 Se3 and Bi 2 Te3. 
The metallic surface states are protected by time-reversal symmetry [3H7] and display the 
linear E{k) dispersion of Dirac fermions as demonstrated by angle-resolved photoemission 
spectroscopy (ARPES) [HI [9] . By the same method, the surface states prove robust [T0HT5] . 
The peculiar spin topology is of helical type and can also be resolved by ARPES when spin 
polarimetry is added to the experiment [To] . 

However, spin-resolved ARPES is a rather demanding detection method since the spin 
polarimetry reduces the ARPES signal by 2 to 3 orders of magnitude. On the other hand, 
there are many interesting subjects that require the measurement of the spin polarization 
and momentum-dependent spin direction. For example, the questions whether the in-plane 
spin component remains locked perpendicularly to the momentum in a warped Dirac cone 
or is in-plane or out-of-plane tilted [T6HT8] . what the maximum value of the polarization is 
[T9H23] . and how the spin behaves in the subsurface region [23] and during hybridization 
with bulk states [2"5] . 

It is, therefore, attractive to search for an alternative method of investigation that gives 
the same information. In photoemission of core levels, the spin resolution is largely com- 
plementary to circular dichroism in which the photoemission intensities for excitation with 
left- and right-circularly polarized radiation, Il(E) and Ir(E), respectively, are evaluated 
as a function of binding energy E. The data from both methods can be analyzed straight- 
forwardly with the same atomic model, as the example of spin-split 5p core levels of Gd and 
Tb shows [261 127] . Also in the valence band, a circular dichroism effect in ARPES is present 
and can be related to the electron spin in ferromagnetic transition metals [2H1 EH] • 

We have recently observed a strong circular dichroism effect in ARPES from the topo- 
logical surface state of Bi 2 Te3 [21J. Spin-resolved ARPES from the same system excited 
by linearly polarized light has been compared to the circular dichroism in ARPES without 
spin resolution. This showed that the circular dichroism reverses together with the spin 
texture at the binding energy of the Dirac point [21] . The circular dichroism asymme- 
try A = (I L — Ir)/(Il + Ir) was found to be very large (> 20% [21]) when compared to 
measurements from ferromagnets (~ 3 to 5% [28, 29J). 

For Bi 2 Se 3 , a strong spin polarization of the topological surface state has been observed 
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by spin-resolved ARPES [18], [22j [23], and in addition several measurements of the circular 
dichroism in ARPES have been reported for this system [3U1 EU [33J . Wang et al. conducted 
circular dichroism measurements at a photon energy of hu = 6 eV using a pulsed laser source 
and a time-of- flight detector [30]. The dichroism effect is interpreted assuming transitions 
into a spin-degenerate continuum of final states and, therefore, is sensitive only to the spin 
in the initial state of the photoemission process. By measuring one spin component in 
the surface plane under two different angles and applying symmetry arguments, the two in- 
plane spin components (S x ), (S y ) and the perpendicular component (S z ) are determined [30J. 
Another measurement on Bi 2 Se3 at photon energies of 10 and 13 eV led to similar results and 
a circular dichroism effect of 30% [31] . Park et al. concluded that for these photon energies a 
free-electron final state can be assumed and for left- and right-circularly polarized light final 
states of different orbital angular momentum character are reached. The orbital angular 
momentum was found to be locked to the momentum in a similar way as the spin, and 
the orbital and spin angular momenta were determined to be antiparallel to each other [3T] . 
Jung et al. reached the same conclusion of antiparallel orbital and spin angular momenta for 
Bi 2 Te3 [32]. Ishida et al. compared Cu-doped Bi 2 Se3 to SrTiC-3 for which a dichroism effect 
> 60% is observed [33]. The measured circular dichroism of Cu x Bi 2 Se3 at hu ~ 7 eV appears 
similar to that reported in the other studies but is assigned to a geometrical origin. In the 
corresponding experimental setup there are two planes perpendicular to the sample surface 
for which the dichroism disappears [33]. This is first of all the detection plane spanned by 
the incident light and the photoelectron momentum. Furthermore, another plane normal 
to the surface and perpendicular to the detection plane is identified where the dichroism 
disappears as well because of the wave function symmetry. This plane is a nodal plane 
only for circular dichroism of the topological surface state due to a mirror symmetry of 
its effective Hamiltonian and its two-dimensionality, in contrast to the bulk valence band 
states. Accordingly, it can be seen as a specific feature reflecting the two-dimensionality of 
the electronic states [33]. 

In the present work, we investigate the circular dichroism in ARPES of Bi 2 Te3 for different 
photon energies and demonstrate that our data question the previous interpretations. It is 
concluded that current explanations are too simple and a more detailed theoretical treatment 
is required. Corresponding calculations are conducted and compared to the experiment. 

We have grown single crystals of Bi 2 Te3 by the Bridgman method and cleaved them in 
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situ. The achieved (111) surfaces are of high quality as concluded from the sharp features in 
angle-resolved photoemission of the valence band (Fig. 1). Measurements have been carried 
out in ultrahigh vacuum of 1 x 1CT 10 mbar at low temperature (30-40 K) with a Scienta R8000 
electron analyzer at the UE112-PGM2a beam line of BESSY II with circularly-polarized 
undulator radiation. The geometry of the experiment is shown in Fig. 1. 

The bulk electronic structure of Bi 2 Te3 is obtained by performing first-principles calcu- 
lations within the framework of the density functional theory using the generalized gradient 
approximation to model exchange-correlation effects [M] . The spin-orbit coupling is included 
in the self-consistent cycles of the electronic structure calculation. 

The angle-resolved photoemission intensity calculations are based on the one-step model 
[55] . We use a fully relativistic formalism, allowing to consider in a natural way effects in the 
photocurrent calculation induced by spin-orbit coupling because the practical calculation is 
based on the Dirac formalism [361 EZ] • The spin-orbit coupling enters thus in the calculation 
of the ground state and again in the photoemission calculation. The photoemission calcu- 
lation itself is based on multiple-scattering theory, using explicitly the low-energy electron 
diffraction (LEED) method to calculate the initial and final states for a semi-infinite atomic 
half-space. In this way the final state is calculated by the best available single-particle 
approach as a so-called time-reversed LEED state [38]. In line with this, the initial state 
is represented by the retarded one-electron Green function for the same semi-infinite half- 
space. The photoemission calculations include matrix-element effects, multiple scattering 
effects in the initial and final states, the effect of the photon momentum vector, and the 
escape depth of the photoelectrons via an imaginary part in the inner potential. These 
lifetime effects in the final states have been included in our analysis in a phenomenological 
way using a parameterized complex inner potential V (E) = V m (E) + iVf(E). Herein, the 
real part serves as a reference energy inside the crystal with respect to the vacuum level. 
To account for impurity scattering, a small constant imaginary value of V\ = 0.004 eV was 
used for the initial state. A realistic description of the surface potential is given through 
a Rundgren-Malmstrom barrier [3H] which connects the asymptotic regime z < za to the 
bulk muffin-tin zero V or by a third order polynomial in z, spanning the range za < z < ze- 
In other words, za defines the point where the polynomial region starts whereas ze defines 
the point where the surface region ends and the bulk region starts with the first atomic 
layer. The effective ^-dependent surface barrier V is scaled with respect to the vacuum level 
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E vac = 0.0 eV utilizing the value of the work function = 5.0 eV. The zero of the z scale 
lies in the uppermost layer of atoms. 

The geometrical origin for a circular dichroism in ARPES as has been described in detail 
in Ref. [29] can be understood on the basis of Fig. 1. No dichroism is expected when 
the geometry is such that a symmetry operation which transforms right- into left-circularly 
polarized light leaves the momentum vector of the photoelectron k unaffected. This is shown 
at the top of Fig. 1. The green plane (77 = 0°, k y = 0) is the plane of incidence and is a mirror 
plane transforming the right circularly polarized light of Fig. 1(a) into the left circularly 
polarized light of Fig. 1(b) without affecting k. The green plane constitutes, therefore, a 
nodal plane meaning zero dichroism for k x = 0. For all other situations a dichroism can 
occur, for example in the yellow plane (k y = 0) which is the detector plane. (The detector 
slit is indicated as a narrow vertical rectangle in Fig. 1.) Reversal of the light polarization 
changes k from +1] to —7] [Fig. 1(c)] which means opposite sign of the dichroism for emission 
angles +r] and —r\ [Fig. 1(d)]. Figure l(e)-(g) show separately the photoemission intensity 
for right- and left-circularly polarized light of Bi2Te3. The effect is apparently very large 
and visible already comparing the intensities 1l and Ir. The corresponding asymmetry A is 
then shown in a color representation in Figs. 2-4. Figures 2-4 display a white line in their 
center, k y = 0, that represent the nodal plane mentioned above. As can be seen, not only the 
topological surface state shows the nodal plane but also the bulk states, a fact which would 
support the interpretation as a geometrical effect. However, it should be stressed that the 
additional perpendicular nodal plane discussed in Ref. [33] is identical to the plane in which 
the data of Fig. l(e)-(g) and Fig. 2 have been measured. This means that the description 
suggested by Ishida et al. as dichroism due to geometrical effects only is not appropriate 
for our present case. Moreover, the reversal of the dichroism effect at the binding energy 
of the Dirac point contradicts the model by Ishida et al. because this behavior cannot be 
explained based on purely geometrical effects. 

Figure 2 shows the behavior of the circular dichroism in the experiment at photon energies 
between 21 and 100 eV. The asymmetry is very large (80% at hv = 55 eV). (Note that, as a 
two-dimensional state, the topological surface state has a binding energy independent of the 
photon energy which may still change slightly within Fig. 2 due to surface doping by residual 
gas.) The top row [Figs. 2(a-d)] shows complete £7(kii) dispersion relations which are all 
characterized by a circular dichroism effect that changes sign with binding energy at the 
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Dirac point. Generally, the lower-binding-energy range shows a large dichroism effect, (ky 
is the projection of the electron wave vector on the surface plane.) This range contributes 
spectral intensity of the topological surface state only while the higher-binding energy range 
has to some extent bulk-like contributions as well. We see that for each selected photon 
energy in Figs. 2(a-d) the complete circular dichroism signal has reversed sign, above as 
well as below the Dirac point. Figure 2(e) shows for hv < 70 eV the detailed behavior with 
smaller steps in photon energy. The behavior around 45 eV may indicate additional sign 
changes [40J. According to the initial-state model for the spin [30], the data in Fig. 2 would 
mean that the spin of the probed ground state changes during the photon energy scan. This 
is not possible, not even under the assumption of a layered spin texture, which has been 
predicted to reverse in the topmost atomic layers [21], because the photocurrent from the 
topmost layer dominates in the photoemission signal from the topological surface state. A 
similar problem arises with the initial-state orbital angular momentum model [31] because, 
again, it is unclear how the probed initial state can depend on the photon energy of the 
probing radiation. 

We have, therefore, performed calculations for the ground state and calculated the pho- 
toemission spectra as described briefly above. Figure 3(a) shows the resulting dichroism in 
ARPES at hv = 27 eV using the same representation as for the experimental data. As can 
be seen, the angle- and the binding-energy-dependences are very similar to the experiment 
in Fig. 2. For investigating the origin of the dichroism effect further, our photoemission 
calculations allow us to vary the spin-orbit-coupling strength. The idea behind this is that 
the decoupling of spin and orbital moment will reduce the effect that the spin polarization 
can have on the circular dichroism. A reduced spin-orbit coupling can be simulated simply 
by an increased speed of light cq. Before applying this in Fig. 3, we tested the use of 1.09co 
and 1.24co on the spin-orbit splitting of W metal between and r§ + states and obtained 
90% and 75%, respectively, of the original splitting. This effect is not too far from the 
expected scaling with 1/cq in the atom. The first impact that one notes in Fig. 3 is the 
change in the binding energy of the Dirac point. This is due to the effect of the spin-orbit 
coupling on the bulk band inversion. It should be stressed that if the spin-orbit coupling 
would be reduced further the band inversion and with it also the topological surface state 
would disappear. Figures 3(b) and 3(c) show that the changes of the circular dichroism of 
the topological surface state are rather small when the spin-orbit coupling is reduced. This 
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suggests that the contribution of the spin to the circular dichroism is a minor one, and we 
will return to this question further below in connection with the photon-energy dependence. 

One-step-photoemission intensities in general are based on ground state electronic struc- 
ture calculations and, as a consequence, the energetics of the final states obtained for higher 
photon energies often deviates from the experimental situation depending on the excitation 
energy. In this sense, our calculations have qualitative character and do not, e. g., reproduce 
the position of the experimental sign change between 21 and 25 eV photon energy. However, 
Figure 4 shows that between 25 and 50 eV the sign has clearly reversed in the calculation. 
This is an important confirmation of the experimental results of Fig. 2 demonstrating that 
the sign change is reproduced by our photoemission theory and can be ascribed to the final 
states. 

We have previously investigated the L-gap surface state of Cu(lll) by one-step photoe- 
mission calculations for a wide photon-energy range from 21 to 70 eV and compared to 
ARPES experiments |4I] . In that case, the dichroism asymmetry depends on the photon 
energy as well. An analysis of the final states confirmed that the dichroism effect is strong 
where a <i-type final state is reached in agreement with the expectations from selection rules 
for the orbital angular momentum [UJ. This mechanism has been discussed later on in a 
very similar context in an ARPES study on Cu(lll) and Au(lll) [42J. The same mecha- 
nism underlies the present results which also involve transitions from p-type initial states 
to d-type final states. When the present calculation is modified to exclude transitions into 
G?-type final states, the sign change between 25 and 50 eV photon energy disappears [40J. 
The final states are, however, more difficult to analyze than in fee Cu due to frequent back- 
folding because of the small size of the bulk Brillouin zone along z (perpendicular to the 
surface). An analysis of the initial state identifies the topological surface state as being due 
to all three p-orbitals, and this is the reason for the stronger dependence on the final states 
concerning the sign changes as compared to the p z -type surface state of Cu(lll). The spin- 
orbit coupling is much weaker in Cu when compared to the present system Bi 2 Te3. The fact 
that spin-orbit coupling is not a precondition for a strong circular dichroism in ARPES has 
been demonstrated also early on: In graphite, transitions from the 7r-band into d final states 
lead to dichroism asymmetries of up to 60% at a negligibly small spin-orbit coupling [43]. It 
should be stressed that spin-orbit interaction also induces a k-dependent spin-polarization 
of the initial states, simply due to the presence of the surface [44] . This effect, which is 
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quantitatively considered in our analysis, has to be accounted for to fully understand the 
experimental data. 

Figure 4 also allows us to inspect the potential of spin-resolved photoemission and its 
depencence on the final states. The bottom row of Fig. 4 shows the calculated spin- 
polarization P = (Jr — 7^) / (JT + 1^) of the photoelectrons as a function of the photon energy 
for linearly polarized light. The sign of the spin polarization remains constant and also its 
value is practically the same (P = 80% at 25 eV and 75% at 50 eV when evaluated 100 meV 
above the Dirac energy). This means that spin-resolved photoemission is much less affected 
by final-state effects than the circular dichroism and, therefore, can be used to deduce the 
electron spin in the initial state rather directly. This has to do with the fact that the circular 
dichroism is the difference of two signals and by itself very sensitive while for the initial-state 
spin the phase space for excitation by linearly polarized light is less restricted. 

In summary, we have investigated the dependence of circular dichroism in ARPES from 
Bi 2 Te 3 on the photon energy. Even for final-state energies that are 10 times higher than 
those considered free-electron like, we observe a sensitive dependence on the photon energy 
and a reversal of the sign. The various proposed initial-state models favoring the spin or 
the orbital angular momentum cannot be applied, and the same holds for pure geometric 
models. The circular dichroism in ARPES has clearly been identified as a final-state effect 
experimentally with confirmation by the results of one-step photoemission calculations. The 
spin polarization of the photoelectrons is not affected. 

O. R. thanks G. Bihlmayer for helpful discussions. Financial support from the Deutsche 
Forschungsgemeinschaft (Grant No. EB-154/18) and the Bundesminsisterium fur Bildung 
und Forschung (Grants No. 05K10WMA) is gratefully acknowledged. 

Note added: Recent one-step photoemission calculations [45] do not support the con- 
clusions by Jung et al. [32]. A normal-incidence geometry is suggested as probe for the 
initial-state spin perpendicular to the surface [45]. 

*Present address: Physikalisches Institut, Universitat Wiirzburg, Am Hubland, 97074 
Wiirzburg, Germany. 
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FIG. 1. (Color online) (a-d) Experimental photoemission geometry and comparison between 
Bi 2 Te 3 data along the T-K direction excited by (e) linearly and (f) right and (g) left circu- 
larly polarized light of hv = 55 eV. The topological surface state is identified by its linear 
dispersion upwards from the Dirac point around 0.25 eV binding energy, and the dichroism 
is strongly visible without need for difference spectra. 

FIG. 2. (Color online) Change of the circular dichroism effect with photon energy, (a-d) 
The upper row shows three reversals of the sign between photon energies of 21 and 100 eV. 
(e) Plot of the dichroism asymmetry A evaluated 100 meV above the Dirac point vs. photon 
energy. The dichroism asymmetry is large (80% at hv = 55 eV). 

FIG. 3. (Color online) Calculations for a reduced spin-orbit interaction for the example of 
hv = 27 eV. The spin-orbit coupling is reduced from (a) 100% (c ) to (b) 90% (1.09c ) 
and (c) 75% (1.24c ). This modifies the inverted bulk band gap and moves the topological 
surface state in energy. The circular dichroism remains very similar indicating a minor role 
of the electron spin polarization for the circular dichroism. The overall change in binding 
energy is a side effect and not relevant here. 

FIG. 4. (Color online) Results from the one-step photoemission calculation. Top: The 
calculated circular dichroism changes sign between photon energies of (a) 25 eV and (b) 
50 eV. Bottom: For the same system but linearly polarized light. The calculated spin 
polarization P of the photoemission from the topological surface state is unaffected by the 
photon energy. It reaches P ~ 80% at 25 eV in (c) and 75% at 50 eV in (d). 
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